Multi-photon, multi-mode polarization entanglement in parametric down-conversion 
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We study the quantum properties of the polarization of the light produced in type II spontaneous 
parametric down- conversion in the framework of a multi-mode model valid in any gain regime. We 
show that the the microscopic polarization entanglement of photon pairs survives in the high gain 
regime (multi-photon regime), in the form of nonclassical correlation of all the Stokes operators 
describing polarization degrees of freedom. 
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I. INTRODUCTION 

The quantum properties of light polarization have been widely studied in the regime of single photon counts. In 
comparison, only recently there has been a rise of interest towards the quantum properties of the polarization of 
macroscopic light beams 0, 0, 13- E| j mainly due to their potential applications to the field of quantum information 
with continuous variables and to the possibility of mapping the quantum state from light to atomic media Q- 
A well-known source of polarization entangled photons is parametric down-conversion in a type II crystal. Here, a 
pump field at high frequency is partially converted into two fields at lower frequency, distinguished by their polar- 
izations. Due to spatial walk-off in the crystal, the two emission cones are slightly displaced one with respect to the 
other, and the the far-field intensity distribution has the shape of two rings, whose centers are displaced along the 
walk-off direction, as e.g. shown by Fig.l. The two regions where the far-field rings intersect have a very special role. 
In the regime of single photon pair detection, the polarization of a photon detected in one of this region is completely 
undetermined. However, once the polarization of one photon has been measured, the polarization of the other photon, 
which propagates at the symmetric position, is exactly determined. In other words, when considering photodetection 
from these regions, the two-photon state can be described as the ideal polarization-entangled state[g- Photons pro- 
duced by this process has become an essential ingredient in many implementations of quantum imformation schemes 
(see e.g. 0,0). 

The question that we address in this paper is whether the microscopic photon polarization entaglement leaves any 
trace in the regime of high parametric down-conversion efficiency, where the number of down-converted photons can 
be rather large 9], and in which form. 

To this end parametric down-conversion is described in the framework of a multi-mode model, valid for any 
gain regime, which includes typical effects present in a realistic crystal, as diffraction and spatio-temporal walk- 
off. Quantum-optical polarization properties of the down-converted light are described within the formalism of Stokes 
operators. These operators obey to angular momentum- like commutation rules, and the associated observables are 
in general non compatible. We define a local version of Stokes operators and study the quantum correlation between 
Stokes operators measured from symmetric portions of the beam cross-section in the far field zone. In the regions 
where the two down-conversion cones intersect we find that all the Stokes operators are correlated at the quantum 
level. Although the light is completely unpolarized and Stokes operators are very noisy, a measurement of a Stokes 
parameter in one of these regions in any polarization basis determines the value of the Stokes parameter in the 
symmetric region within an uncertainty much below the standard quantum limit. 

A continuous variable polarization entanglement, in the form of quantum correlation between Stokes operators of 
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FIG. 1: Parametric down-conversion from a type II crystal showing the two down-conversion cones at degeneracy. 
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two light beams, have been recently demonstrated 3]. In this work the entanglement is of macroscopic nature, and 

spatial degrees of freedom do not play any role since the beams are single-mode. Continuous variable polarization 

entanglement which takes into account spatial spatial degrees of freedom of light beams is described in [Tj| , where we 

study the properties of the light emitted by a type II optical parametric oscillator below threshold. 

The analysis of this paper is rather focussed on providing a bridge between the miscroscopic and macroscopic domain, 

since our model is able to describe the polarization entanglement in parametric down-conversion with a continuous 

passage from the single photon pair production regime to the regime of high down-conversion efficiency 

Besides its fundamental interest, we believe that the form of entanglement described in this work can be quite promising 

for new quantum information schemes, due the increased number of degrees of freedom in play (photon number, 

polarization, frequency and spatial degrees of freedom), and is well inserted in the recent trend toward entangled state 

of increasing complexity (see e.g. [Toj . where a four-photon polarization entangled state is characterized) 

The paper is organised as follow. Section [H] describes the model for spontaneous parametric down-conversion, in 
terms of propagation equations for field operators. Similar models are known in literature (see e.g. |13| and references 
quoted therein), but besides presenting it in a systematic way, we include all the relevant features of propagation 
through a nonlinear crystal and we provide a precise link with the empirical parameters of real crystals. Section IlIII 
is devoted to the description of the quantum polarization properties of the down-converted light. Stokes operators 
definition and properties are briefly reviewed in section IIII Al In Sec IIII Bl we generalize this definition to a local 
measurement in the beam far field plane, and we introduce the spatial correlation functions of interest. Analytical 
and numerical results for the degree of correlation of the various Stokes parameters detected from symmetric portions 
of the beam cross-section are presented in sections IIII CI IIII Dl both in the case when a narrow frequency filter is 
employed fSec llll iTTIl . and when the filter is broad-band fSec lllllT2"|) . Section HVI provides an alternative description 
of the system and of its polarization correlations in the framework of the quantum state formalism. Section Ivl finally 
concludes. 

II. A MULTI MODE MODEL FOR TYPE II PARAMETRIC DOWN-CONVERSION 

A. Field propagation 

The starting point of our analysis is an equation describing the propagation of the three waves (signal, idler and 
pump) inside a nonlinear x' 2 ' crystal. We consider a crystal slab of length l c , ideally infinite in the transverse directions, 
cut for type II quasi-collinear phase-matching. In the framework of the slowly varying envelope approximation the 
electric field operator associated to the three waves is described by means of three quasi-monochromatic wave-packets. 
We take the z axis as the laser pump mean propagation direction (Fig.l), and indicate with x = (x,y) the position 

coordinates in a generic transverse plane. E^ + \z,x,t), i — o,e,p designate the positive frequency part of the field 
operator (with dimensions of a photon annihilation operator) associated to the ordinary (i = o, the "signal") and 
extraordinary (i = e, the "idler") polarization components of the downconverted beam, and (i — p the "pump") the 
high frequency laser beam activating the down-conversion process. Next we introduce their Fourier transform in time 
and in the transverse domain: 

Mz,q,Q) = lf-l -^=e-^ s e^+^El + \z,x,t) i = o,e,p (1) 

J 2tt J V27T 

Here q is the transverse component of the wave- vector and f2 represents the frequency offset from the carriers uj j ruj e = 
uj p . In the following, we shall assume degenerate phase matching, so that lo — u e = uj p /2. It is convenient to subtract 
from the field operators the fast variation along z arising from linear propagation inside the birefringent crystal. We 
write: 

Ai(z, q, 0) = exp [ik lz (q, S})z]ai(z, q, Q) , (2) 

where k iz (q,Q) — U (q, wT+Oj — q 2 is the projection of the wave- vector along the z direction, with ki(q,u>i + il) 
being the wave number of the i-th wave. In the absence of any nonlinear interaction, we would have 

^-ai(z,q,Sl) = (3) 
az 

being Eq.Q with a,i(z,q, Q) — h{{z — 0, q, f2) the forward solution of Maxwell wave equation in linear dispersive 
media. For the pump wave, we assume that the intense laser pulse is undepleted by the down-convertion process, so 
that a p (z, q, SI) = a p (z = 0, q, Q). Moreover, we assume that the pump is an intense coherent beam and the operator 
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can be replaced by its classical mean value a p {q, SI). 

For the signal and idler beams, the variation of di operators along z is only due to the nonlinear term, proportional 
to the x material second order susceptibility. This is usually very small, so that aj are slowly varying along z. 
This allows us to neglect the second order derivative with respect to z in the wave-equation. Hence the resulting 
propagation equation takes the form (see also for more details, and 0] f° r an alternative derivation): 

(z, q, St) = X J dq' J dSl' a p (q + q' , SI + Si') a] (z, q', Si') e -^n^M,/l c i^j = ,e, (4) 

where % is a parameter proportional to the second order susceptibility of the medium, and 

A %3 (q, q"; SI, Si') = l c [k tz (q, SI) + k jz (<?', SI 1 ) - k pz (q + q" , SI + SI')} (5) 

is the phase mismatch function. Equation 0) describes all the possible microscopic processes through which a pump 
photon of frequency lj p + SI + Si' , propagating in the direction q + q' is annihilated at position z inside the crystal, 
and gives rise to a signal and an idler photon, with frequencies uj p /2 + SI, uj p /2 + Si', and transverse wave vectors q, 
q', with an overall conservation of energy and transverse momentum. The effectiveness of each process is weighted by 
the phase mismatch function J^J, which accounts for conservation of the longitudinal momentum. In the limit of an 
infinitely long crystal, where longitudinal radiation momentum has to be conserved, only those processes for which 
Ajj = are allowed. For a finite crystal, however, the phase matching function has finite bandwidths, say qo in the 
transverse domain and SIq in the frequency domain. 

Equation |@j ) couples all the signal and idler spatial and temporal frequencies within the angular bandwith of the 
pump 5q w -j-, with w p being the pump beam waist, and within the pump temporal spectrum SSI « 1/t p , where t p 
is the pump pulse duration. In general, no analytical solution is available and one has to resort to numerical methods 
in order to calculate the quantities of interest, as described injllj. 

A limit where analytical results can be obtained is that of a pump waist and a pump duration large enough, so that 
5q << qo, SSI << SIq- In this case the pump beam can be approximated by a plane wave 



Equation reduces to 



a p (q + q*,Sl + SI') -> a p S (q + q')5(Sl + SI') , (6) 



l c ^-a {z,q,Sl) = aal(z,-q,-Sl)e- iA ^' n >^ 



d 



z 



l c ^a e (z, -q, -SI) = aat (z, q, SI) e - iA ^^ , (7) 

where a = l c X a p is a linear gain parameter, and 

A(£ SI) = l c [k oz {q, SI) + k ez {-q, -SI) - hp] • (8) 

is the phase mismatch of a couple of ordinary and extraordinary waves propagating with symmetric transverse wave 
vectors q and —q, and with frequencies tv p /2 + SI, lu p /2 — Si. 

Solution of the propagation equation (J7J is found in terms of the field distributions at the input face of the crystal. 
Coming back to the field operators defined by Eq. 0), we define the field operators at the input and output faces of 
the crystal slab as 

AT(q,Sl) = ai(z = 0,q,Sl) (9) 
A° ut (q,Sl) = ai(z = l c ,q,Sl)exp{ik iz (q,n)lc} (10) 

By solving Eq.lQ) the transformation from the input to the output operators is found in the form of a two- mode 
squeezing transformation: 

A™\q,Sl) = U (q,n)A™(q,n) + V (q,n)At™(-q,-Sl) 

A° ut (q, n) = U e (q, Sl)AT(q', SI) + VM V)At in (-q, -SI) , (11) 

linking only symmetric modes q, SI and —q, —SI in the signal and idler beams (see e.g. [T^ | for a similar transformation 
in the type I case) . If we require that free space commutation relations 

A?(q, SI), Af{q', SI')} = 5 id 5(q- q')6(Sl - Si') i,j = o, e (12) 
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are preserved from the input to the output, it can be easily shown that the complex coefficients of the transformation 
need to satisfy the following conditions: 

IW,fi)| 2 -|t$(?,fi)| 2 = 1 (i = o,e) (13) 

u {$,n)v e (-q,-n) = v (<f,n)u e (-?,-(i) (u) 

By taking the modulus of the second relation and making use of the first two ones, the complex equation (|14(l can be 
written as two equivalent real equations: 

\V (q,il)\ 2 = \V e (-q,-Q)\ 2 (15) 

arg [U Q {q, 9)V e {~q, -O)] - arg [V (q, Sl)U e {-q, -fi)] := 2^(q, fi) (16) 

With this in mind, the coefficients of the transformation llljl can be recasted in the form 



U (q, fi) = U(q, Siy*^ , V (q, SI) = V(q, n)^®") 

U e (q, n) = U(-q, -0) e - i * , (-«'- n ) , V e (q, n) = V(-£ -0)e- i v(-9*- n ) 



(17) 



with 



f/(£Q) = coshr(g,fi)e^ ( «< n) e ie( «< 0) 

V(£fi) - sinhr(g,fi)e^® n )e- ie ^ n ), (18) 

where r(<f, Q),(p(q, f2, ^(g, f2, 6(q, f2, are independent real functions of q, fl. 

We outline that the form of the transformation ((11(1 , together with the unitarity requirements ((131 114(1 , are enough 
the derive the general form of the results presented in this paper. In the following, we shall present results for a 
specific device, namely travelling wave parametric down-conversion, and we shall take as an example the case of a 
BBO (6eta-barium-borate) crystal. However, a similar investigation can be carried out for any device characterized 
by an input/ouput transformation of the form 1(111) . The case of Type II parametric down-conversion inside an optical 
resonator is for example investigated in [l4j . 

More insight into the problem is gained by looking at the explicit solution of the propagation equation (JTJ. We 
obtain 

U(q, Q) = jcosh [T(q, O)] + i^pjy sinh [T(q, «)]} , (19) 

V(q,Q) = e^-^— sinh[r(q,fi)] , (20) 
T(q,n) 

<p(q, SI) = | [KM fi) - k ez (-q, -fi)] , (21) 

with 



r®n) = do>-*4W, (22) 



and A(tf, fi) is the phase mismatch function defined by 



B. Phase matching curves 

The gain functions 1(191 120(1 reach their maximum value for phase matched modes, that is, the modes for which 
A(<f, fi) = 0. By assuming the validity of the paraxial and slowly varying envelope approximations, the longitudinal 
wave- vector components ki Z (q, fi) can be expanded in power series of <f, fi. By keeping only the leading terms we 
obtain: 

1 ld 2 fc dk- a 2 

MM « *. + + ^ 2 + d£* " & 4 = °' e ■ < 23 > 

The first term at r.h.s. is fcj = Jij(k>j, <f = 0) cjj/c, with being the index of refraction at the carrier frequency of an 
ordinary (extraordinary) wave propagating along z direction. The second term, ~fl = ^jf-fi, accounts for the fact 
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that the three wave-packets move with different group velocities v % . The third term describes the effects of temporal 
dispersion. In writing the fourth term, we assumed that the crystal is uniaxial and the crystal optical axis lies in the 
z-y plane. This term is present only for the extraordinary waves , and ^p 1 = — pi where pt is the walk-off angle of the 
wave. Finally, the last term describes the effects of diffraction for a paraxial wave. 
With this in mind, the phase matching function can be written in the form: 



a 2 I 

A(g, fi) = A + p e l c q y - + flr coh + -e{VLT coh ) 2 (24) 

9o 1 



where 
1) 



A = (k + k e - k p )l c (25) 

is the collinear phase mismatch (i.e. the phase mismatch of the three waves at the carrier frequencies when 
propagating along the longitudinal direction); 



2) 



3) 



^ _ / 1 k e + k _ 2tt n e + n a 
V l c 2k e k Q V A/ c 2n e n 

with A = Anc/ujp being the wavelength in vacuum at the carrier frequency lo p /2, and n e ,n a the odinary and 
extraordinary refraction indexes inside the crystal at the carrier frequency. This parameter defines the typical 
bandwidth of phase matching in the transverse q-space domain. Its inverse l co h = 1/qo will be referred to as 
the coherence length. 



r coh = l - o - l - e (27) 

v° v% 

with v l g being the group velocities of the two waves, is the the difference between the time taken by the signal 
and idler wave-packets to cross the crystal. This defines the typical scale of variation of gain functions in the 
temporal domain for type II phase matching, and it will be referred to as the amplifier coherence time. 



4) Finally 



/d 2 fc d 2 k e \ l c 
{-dW + dw) r^ oh 



e =(^ + ^)r^ (28) 



is a dimensioneless parameter that depends on the temporal dispersion properties of the signal and idler pulses 
(typically e << 1). 

The equation A(<f, O) = defines in the (q x , q y ) plane a circonference, centered at the position 

q x = 0, q y = qc = T^qlPelc (29) 

and with radius given by 

I ~q 2 1 
qn = qo\ A + — + flTcoh + -ze{VlT coh ) 2 . (30) 
V % * 

This corresponds to the phase matched modes for the signal (ordinary) wave. Phase matched modes for the idler 
wave, emitted at the frequency — f2, lye on the symmetric circonference. 

Figure plots some examples of this phase matching circles, in the form of polar plot, with 9 being the polar angle 
from the pump direction (z axis) outside the crystal and 4> the azymutal angle around z. Parameters are those of 
a BBO crystal, cut for degenerate phase matching at 49.6 degrees, for a pump wavelength of 351 ran. They have 
been calculated with the help of empirical Sellmeier formulas for refraction indices in Ref. 15|- For comparison, 
superimposed to the curves calculated b y m eans of Ea. (|24|) . the figure shows the "exact" phase matching curves, 
calculated with the method described in jl6| | , by means of a public domain numerical routine available at [17| . The 
plots show a rather good agreement, in any case within the error inplicit in the use of empirical Sellmeier formulas. 
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signal 702 nm 




FIG. 2: Polar plot of phase matching curves in a BBO crystal. Comparison of the approximated formula 1241 with the "exact" 
phase matching formula, calculated with the method in Hal (NIST). 6 is the polar angle from the pump direction of propagation, 
cj> is the azymuthal angle. X S i g nai = tidier = 702nm. (A pump = 351nm), and the pump propagates at an angle of 49.6 degrees 
from the crystal optical axis 



III. POLARIZATION CORRELATION: QUANTUM FIELD FORMALISM 



A. Stokes operators:definition and properties 



Quantum-optical polarization properties of light are conveniently described within the formalism of Stokes operators, 
which represent the quantum conterparts of the Stokes vectors of classical optics. The polarization state of a classical 
beam can be described by means of a Stokes vector, and of its associated Poincare sphere. Stokes vectors pointing on 
the equator of the sphere represents linearly polarized light; if the vector points in the positive (negative) directon of 
Si the light is horizontally (vertically) polarized, while the S 2 direction identifies light polarized at 45 (-45) degrees. 
The 5*3 direction corresponds to circulary (right and left) polarized light. A fourth parameter, Sq is the total beam 
intensity, and gives the radius of the Poincare sphere. For a polarized beam +S% + S 3 = Sg, so that the polarization 
state is represented by a point on the sphere surface. 

In the quantum mechanical description of light polarization, Stokes parameters are replaced by a set of four Stokes 
operators. For a single mode of an electromagnetic field, they are defined in terms of the photon annihilation operators 
for a vertical and horizontal linear polarization mode an, «y, as: 

So = a H a,H + ayav (31) 

S\ — a) H a,H — OyCLy (32) 

5*2 = a) H a v + a) v a H = a\ 5 a 45 - a^ 45 a_ 45 (33) 

S 3 = -i (a^ H a v - o^clh) = a* R a R - a[a L , (34) 

where S45 , a_45 denote annihilation operators on the oblique polarization basis, and an are annihilation operators 
on the circular right and circular left polarization basis. The first two operators represent, respectively, the sum and 
difference of photon numbers in the vertical/horizontal polarization basis. Operators S 2 and S3 are the difference of 
photon numbers in the oblique and circular polarization basis, respectively. All these observables can be measured 
by means of a polarizing beam splitter and quarter and half wave plates, as e.g. described in |Tsj . However, while 
operator So commutes with all the others, the remaining three do not: 

[S 1 ,S 2 ]=2iS 3 , [S 2 ,S 3 ]=2iS 1 , [S 3 ,S 1 }=2iS 2 . (35) 

The set of Stokes operators has angular momentum-like commutation relation, and the associated observables are in 
general non compatible. The quantum state of a light beam cannot be any more visualized as a point on the Poincare 
sphere, since quantum noise introduces a minimum uncertainty in the values of the Stokes parameter. Polarization 
squeezed states, whose uncertainty can be represented by an ellipsoid (see e.g. 2}) has been recently realized 0. 
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B. Stokes operator correlation in the far field of parametric down-conversion 



The main idea of this paper is to study the quantum correlation between Stokes operators measured from symmetric 
portions of the far held beam cross-section. To this end, we consider a measurement of the Stokes operators over a 
small region D(x) centered around a position x in the far-field plane of the down-converted field, and over a detection 
time T (tpically we will take T much larger than the crystal coherence time). 

(36) 



where 



Si(x) = dt' dx'a % (x',t') , 

JT JD(x) 

a (x,t) = Al(x,t)A e (x,t) + Al(x,t)A e (x,t) , 
ai{x,t) = Al(x,t)A (x,t) - Al(x,t)A e (x,t) , 
a 2 (x,t) = Al(x,t)A e (x,t) + Al(x,t)A Q (x,t) , 

a 3 (x,t) = -i\Al(x,t)A e (x,t) - Al(x,t)A e (x,t) 



(37) 
(38) 
(39) 
(40) 



A / e denotes the field operator for the ordinary/extraordinary polarized beam in the far-field plane, which can be 

observed in the focal plane of a lens, placed as shown in figure |3J 

By using the free field commutation relations (|12(l , it can be easily shown that 



[Si(x), S 2 (x)} = 2iS 3 (x) , [S 2 (x),S 3 (x)} = 2iS 1 (£) , [S 3 (x), &(£)] = 2iS 2 (x) . 

while operators measured from different (and not connected) detection pixels commute. 
In the following, we shall consider Stokes operator correlation functions of the form: 

(6Si(x) 6Sj(x')) = (Si(x) Sj(x')) - (Si(x)) ( S 3 (x")) , = 0, . . .3) . 

A useful tool for calculation are the correlation functions of the Stokes operator densities H37I40JI 

Gij(x,x';r) = (&i(x, t + r)aj(x' - {&i(x,t + r)) (&j(x',t)) 

and their spectral densities 

G ij (x,x';Q)= [ dTe iQT G lj {x,x';T) 



Their relation with the correlation functions of the measured Stokes operator H42|l is given by: 



{SSii^SSjix')) 



dxi 



D(x) 



D(S') 



T 2 ( T N 
dx[ J dil— sine 2 I Q— ) Gij (xi,x[;Q) 



We notice that: 



T 2 / T 
lim — sine" ( fl 



(41) 

(42) 
(43) 

(44) 
(45) 
(46) 



and that this function acts under the integral as a frequency filter with bandwidth Af2 = 2n/T. We shall assume in 
the following that the detection time is much larger than the crystal coherence time. Under this assumption 



{SSi(x)6Sj(x')) 



T»T coh 



T 



dx\ 



D(x) 



D(x' 



dx[Gij (xi,x[; fl = 0) 



(47) 



When a lens of focal length / is placed at a focal distance both from the crystal output plane and the observation 
plane (se e fig ure ISt, the field operators in the far-field plane are connected to those at the crystal output by the usual 
mapping 20] 



(48) 
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BBO 1 , BBO 2 

V2 @45 



FAR 
FIELD 



FIG. 3: Schematic set-up for a measurement in the far field plane with a compensation crystal. BBOl, down-conversion crystal 
of length l c ; A/2@45 half-wave plate, rotates polarization by 90 degrees; BB02, compensation crystal of length l' c , L lens of 
focal length /. 



where / is the focal length of the lens used to image the far field plane and Ai is the wavelength (in vacuum) at the 
frequency uj p /2 + Sl. 

Since the light statistics is Gaussian (the output operators are obtained by a linear transformation acting on input 
vacuum field operators) all expectation values and correlation functions of interest can be calculated by making use 
of the second order moments of field operators. These can be easily calculated in the far field plane by assuming that 
the down-converted field operators at the input crystal face are in the vacuum state and by using the input/output 
relations toghether with Q48[l . thus obtaining 

^(f, 0)^(5', n')) = Si,j5(x-x')8(n-s}')\Vi(x,n)\ 2 , (49) 

(At&fyAjtf^Sl')) = -^-(l-S hj )S{x^-+x' )S(Sl + Sl')U i (x,Sl)V :i (-x^-,-Sl) , (i,j = o,e). (50) 

\ / Ai ' Ai Ai 

In this formula 

Ui(x, Si) = Ui ($= x^j, Slj , Vi(x, Si) = Vi (q = > *H ( 51 ) 

where U, Vi are are the gain functions defined by (jl7H22(l . It can be noticed the presence of the nonzero "anomalous" 
propagator H50[l , a term which is characteristic of processes where particles are created in pairs. In order to simplify 
the notation, in the following we shall consider the case SI << uj p /2, and take A = A e = A = 2X p . In a real 
experimental implementation, however, the validity of such an approximation should be carefully checked when not 
using narrow frequency filters; twin photons produced at different wavelengths A e , A G , and travelling with symmetric 
q, —q transverse wave vectors are actually propagating at different angles from the pump and will be intercepted in 
the far field at two sligtly different radial positions. 

The fact that the field spatial correlation are perfectly localized in the far field (the Dirac-delta form of the correlation 
peak) is a consequence of the traslational symmetry of the model in the transverse plane (plane wave pump and 
a crystal slab infinite in the transverse direction). A trivial formal fault is that the far field mean intensity of the 
downconverted beams diverges, as a consequence of the infinite energy of a plane- wave pump. This artificial divergence 
can be formally eliminated with the trick used in |l9l |22T | , where a finite size pupil was inserted at the output face 
of the crystal. The spatial Dirac-delta functions in Eas. H49l50|) are substituted by a finite version, and a typical 
resolution area, proportional to the diffraction spot of the pupil in the far field plane, is introduced in the scheme. For 
a pupil of transverse area Sp, this is given by Dp — (Xf) 2 /Sp. The typical scale of variation of the gain functions 
(|51|l in the far field plane is 

X a = q Xf/(2iT) ; (52) 

when Xq is much larger than the resolution area (or, equivalently, when the pupil size is much larger than the amplifier 
coherence length), the mean photon number distribution in the far field plane is given by 



(53) 



(Ni(x)) = / dx' / dt(Al(x" ,t)Ai(x" ,t)) = 

JD(S) JT 

U R JD(x) J W 
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When the finite size of the pump is taken into account in a numerical model |TT| , it is easily seen that the resolution 
area is rather given in terms of the spot size of the pump as it is imaged in the far field plane. For a Gaussian pump 
of waist w p , D R w (A/) 2 / '(ra^). 

In this limit of small resolution area, the mean value of Stokes operators is given by: 



(Si(x)) = 



T 

Dr~ 

T 



dx" 



D{2) 



Dr 



dx" 



R JD(x) 



dn r 

2^ 



\v (x',n)\ 2 + \v e (x',n)\ 2 



\v (x> ,n)\ 2 -\v e (x" ,n)\ 2 



(&(*)) = (S 3 (x)) = o 



(54) 

(55) 
(56) 



C. Correlation in Stokes operators Si, So 



The first and second Stokes operators represent the sum and the difference, respectively, between the number of 
ordinary and extraordinary photons (say horizontally and vertically polarized photons) measured from a detection 
pixel in the far field plane. 

S (x) = N (x)+N e (x) (57) 
Sxix) = N (x)-N e (x) (58) 

The plane wave pump model predicts that the number of ordinary and extraordinary photons collected from any 
two symmetric portions of the far field plane are perfectly correlated observables [ll], [jjj . This result is a direct 
consequence of pairwise emission of photons with horizontal (ordinary) and vertical (extraordinary) polarizations, 
propagating in symmetric directions, as required by transverse light momentum conservation. Hence, this model 
predicts an ideally perfect correlation, both between §o(x), So(—x), and between Si(x), —Si{—x) for any choice of 
the position x in the far field (notice that Sq{x) commutes with S\{5?)). 

In a more sofisticated numerical model |Tlj . it is readily seen that the finite width of the pump profile introduces 
an uncertainty in the directions of propagation of the down-converted photons. As described by the propagation 
equation Q}, when a o photon is emitted in direction q, its twin e photon is emitted in the direction — q within an 
uncertainty Sq cx 2/w p , which is the bandwidth of the pump spatial Fourier transform. A photon number correlation 
well beyond the shot noise level is recovered when photons are collected from regions larger than a resolution area 

In the limit of a small resolution area, long but straigthforward calculations [2l| show that: 

G 00 (x,x';n) = —[6(x-x')Fi(x,n) + S(x + x')F 2 (x,n)] (59) 



Dr 
1 

Dr~ 



Gii (£,£'; ft) = —[6(x-x')Fi(x,n)-S(x + x')F 2 {S,n)] (60) 



with 

Fi(x,Q) = J ^{\V (x,Lu)U (x,Lu + n)\ 2 + \V e (x 1 u;)U e (x,Lu + n)\ 2 } (61) 

F 2 (x,n) = J ^{u {x,Lu)u;(x,Lu-n)v e (-x,-Lu)v:(--x,-Lu + {i) 

+ U e {S,uj)Ut{x,u)-Q)Vo{-S,-u))Va{-S,-uj + n)} (62) 

The correlation functions have two peaks; the first one, located at x' = x, accounts for the noise in the measurement 
of Stokes parameter from a single pixel. The second one is located at x' = —x and accounts for correlation (anti- 
correlation) between measurements performed over symmetric pixels. By taking into account the unitarity relations 
(|13I15|I . it can be immediately noticed that when ft — (corresponding to long detection times) Fi(x, 0) = F 2 (x, 0), 
and the two corelation function peaks have the same size. This represents the maximum amount of correlation allowed 
by Schwarz inequality, which requires that 

i 

(SSi^SSii-x)) < (SSiffiSSifflWSii-aOSSii-x)) * (63) 
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In our case, assuming two symmetric detection pixels D(x) and D{—x), we have e.g 

T 



T 



D, 



R JD(x) 



dx'F^x',0) 



(6Si(g)5Si(-g)) = — ^- / dx'F 2 {x',0) = -(SS 1 (x)SS 1 {x)) 
U R Jd(S) 



(64) 
(65) 



Finally, the existence of such a perfect correlation implies that both Si(x) + Si(—x) and Sq(x) — Sq(—x) are noiseless 
observables. For example: 



( SSi(x) + SSi(-x) ) = 2 (SSi(x)SSi(S)) + (SSi{x)SSi(-x)) 







(66) 



D. Correlation in Stokes operators £2, S3 



Quite different is the situation for the other two Stokes operators S 2 , S3, which involve measurements of the photon 
number in a polarization basis different from the ordinary and extraordinary ones of the crystal, namely in the oblique 
and circular polarization basis. 

Calculations along the same lines of those performed for the first two Stokes operators show that also in this case the 
correlation functions display two peaks, one representing the noise associated to the measurement over a single pixel, 
the other the correlation between symmetric pixels. 



G 22 {x, x';Cl) = G 33 (f,x';Sl) 

= -j- [S(x -x')Hi (2, n) + S(x + x')H 2 (x , Q)] 



with 



ffi(£,fi) 



{ \V (X, Uj)U e (x, OJ + n f + \V e (x, Uj)U {x, UJ + fl)\ 2 } 



H 2 (x,Q) = J —{UZ{£,up B (2,u + (l)VZ{-2,-u)V {-2,-u-Sl) 
+ U*(x,u))U (x,lu + Q)V*(-x, -u)V e (-x, -uj - fi)} 



(67) 
(68) 

(69) 
(70) 



However, at difference with the previous case, the two peaks in general do not have the same size, even for a long 
measurement time. Letting SI = in Eqs. (|69l I70|) and using the definition (|17fl . which is a consequence of unitarity, 
we have 



ffiOM) - j ^{\V(x,oj)U(-x,-uj)\ 2 + \V(-x,-uj)U(x,u J )\ 2 } 

/ J 
-^2Re {U*(x , cj)U(-x, -u)V*{x, u)V(-x, -uj)} , 

Hi(x,0) -H 2 {x,0) = [ — \V*{x,u)U{-x,-uj)- V*{-x,-uj)U(x,uj)\ 2 

. 27T 



where U,V appearing in these equations are the functions defined by Eas. l|19l2t)|l . calculated at q 
Moreover, 

(\s 2 (x) - S 2 (-x)Y) = ( \s 3 (x) -S 3 (-x)] 2 ) 



2T f 

— dx' [H x (x' ,0) - H 2 (x' ,0)} 

U R JD(x) 



(71) 
(72) 
(73) 
x2n/(\f). 

(74) 
(75) 



The noise in the difference between Stokes operators measured from symmetric pixels in general does not vanish, due 
to the lack of symmetry x, ft — * —x, —il in the gain functions. In turns, this reflects the effect of spatial walk-off 
between the ordinary/extraordinary beams (described by the term proportional to q y in the phase mismatch function 
12411 and the group velocity mismatch beteween the two waves (described by the term proportional to fl in I24fl . 
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FIG. 4: a) Far field photon number distribution of the down-converted field. b)Distribution of the noise in the difference 
between §2 measured from symmetric portions of the beam cross-section, scaled to the shot-noise level. The distribution for 
S3 is identical. A step function frequency filter AA = 8 nm wide, centered around the degenerate frequency is used in both 
plots, a = 2 



1. Narrow-band frequency filtering results 
Part (b) of FigU| shows a typical result for the noise in the difference between Stokes operators measured from 

r * 1 2 

small (but larger than Dr) symmetric portions of the far-field. Precisely, the figure shows ( §2(x) — S2(—x) ) = 

1 2 r - 1 

S3 (a?) — S3 (—x) ), scaled to the shot noise level, represented by ( Sq(x) + Sq(—x) ), as a function of the transverse 

coordinate x — (x, y) scaled to Xq. Parameters in this plot are those of a 2 mm long BBO crystal, cut at 49.6 degrees 
for degenerate type II phase matching at 702nm. For comparison, part (a) of the figure shows the mean photon 
number distribution in the far field. The numbers associated with the scale in (a) represent the number of photons 
detected over a resolution area Dr and over a crystal coherence time, that is the mean photon number per mode. Both 
plots have been obtained by filtering the emitted frequencies in a bandwidth AA = 8 nm wide around the degenerate 
frequency, by means of a step function filter. Precisely, we let 



At (x, Q) vWUi (x , fi) + y/1 - f{tt)vi {x, O) % = o, e , (76) 

where Vi(x, Q) are vacuum field operators uncorrelated from the signal and idler fields Ai(x, CI), and the filter function 
f(Q) is this case the step function /(Q) = 1 for Qe[-4p, 4p], f(fl) = elsewhere. 

In plot (b) we see clearly two large dark zones, in correspondence of the intersections of the emission cones, where 
the Stokes operator correlation is almost perfect. Out of these regions, basically no spatial correlation at the quantum 
level exists for Stokes operator S2 and S3. 

Remarkably, at the intersection of the two degenerate emission cones, the light is completly unpolarized. Figure 
shows the distribution of (Si(x)} / (Sq(x)) in the transverse far-field plane, showing that it vanishes at the emission 
ring intersection; recalling that the mean value of S2 and S3 is zero everywhere, this means a vanishing degree of 
polarization in these regions. Moreover, in these regions a measurement of Stokes parameters over a single detection 
pixel is very noisy, as shown by figureEl which plots the distribution of ([SS2(x)] 2 ) = ([(^(x)] 2 ), scaled to the shot 
noise level (Sq(x)}. In this plot the uniform dark background correponds to the shot noise level, while the bright spots 
correspond to a noise level 10-12 times larger than the shot noise. 
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Similar results are obtained in any gain regime. In the small gain limit, the noise statistics associated to a measure- 
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FIG. 5: Degree of polarization of the light downconverted by a BBO crystal. Far field distribution of (Si(x)} / (So(x)) . Same 
parameters as in FigHTI 
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FIG. 6: Noise in the measurement of Stokes parameters of the downconverted light by a BBO crystal. Far field distribution of 
([8S 2 {x)] 2 )/{S (x)} = ([5S 3 {x)] 2 }/{So{x)}. Same parameters as in Figg] 

ment over a single pixel becomes essentially Poissonian, but the correlation between Stokes parameters measured from 
symmetric pixels is basically the same as in the high gain regime. Fig. |7|compares the noise in the difference between 
Stokes parameters measured from symmetric pixels in the small and high gain regimes, plotted as a function of the 
vertical coordinate along the circle of maximum gain for the degenerate frequency. The dashed lines were obtained 
with a = 0.01, corresponding to a mean photon number per mode « 10 -4 , the solid lines with a — 2, corresponding to 
a mean photon number per mode « 15 (see Fig0Ji),. In this plot, the two dark lines are, as usual, obtained by filtering 
the frequencies with a step function (AA = 5nm). For comparison, the two light lines show the results obtained by 
means of a more realistic frequency filter, with a Gaussian profile. Precisely, we take the filter function in Ea. H76|) 
as f(£l) = exp [— (£! 2 41n2)/(AS7 2 )] where Af2 is the full width at half maximum (FWHM), and corresponds to an 
interval in wavelengths of 5nm. In this case, losses introduced by the Gaussian shape of the filter slightly deteriorates 
the correlation. 

The results described above were obtained by exploiting a trick commonly used in the experiments performed in the 
single photon pair regime ( for example in the experiment of |||), in order to partially compensate for the temporal 
and spatial walk-off of the down-converted beams. In the regime of single photon pair production, the ordinary and 
extraordinary photons can be in principle distinguished because of their different group velocities inside the crystal, 
and because of their offsets in propagation directions due to walk-off effects inside the crystal. The mere existence of 
this possibility is detrimental for the entanglement of the state. As it is discussed in details in Appendix A, in the 
general case (arbitrary number of down-converted photons), the group velocity mismatch and the spatial walk-off arc 
responsible for the appearance of a propagation phase factor that lower the value of the correlation function between 
Stokes operators measured from symmetric regions. In principle, this problem can be solved by using a vey narrow 
frequency filter, and by performing the measurement over narrow regions centered around the ring intersections. 
However, this lower the efficiency of the set up. An other possibility is to insert a second crystal, after the pump 
beam has been removed, and after the field polarization has been rotated by 90° (see Fig|3J). In this way, the slow 
and fast wave in the first crystal become the fast and slow wave, respectively, in the second crystal, and the direction 
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FIG. 7: Noise in the difference between 5*2 (5*3) measured from symmetric portions of the beam cross-section, scaled to the 
shot-noise level, as a function of the vertical coordinate y along the maximum gain circle for the degenerate frequency. Dashed 
lines: a = 0.01, solid lines a = 2 Light lines: Gaussian frequency filter, of FWHM = 5nm. Dark lines: step function filter 5nm 
wide. 



of walk-off is reversed. At difference from the single photon pair regime, the correlation is optimized when the length 
of this second crystal is chosen as 

l' c = l C ^, (77) 

where a is the linear gain parameter, proportional to the pump amplitude and to the first crystal length (see Appendix 
A) . The fact that the optimal length of the compensation crystal decreases with increasing gain can be understood 
as following |2.'il |24| : in the regime of single photon pair production (limit a — ► 0), the photon pair can be produced 
at any point along the crystal length with uniform probability, so that the average temporal delay of the two photons 
due to the group velocity mismatch are those corresponding to half of the crystal length, and best compensation is 
achieved for l' c = k^. In the large gain regime, more and more photon pairs are produced towards the end of the crystal 
(the number of down-converted photons increases exponentially with the crystal length) , so that walk-off effects are 
best compensated by a shorter crystal, whose length is given by formula J77J. 
When this kind of optimization is not possible, our calculations show that similar results can be obtained by a 




FIG. 8: Effect of the compensation crystal. Noise in the difference between §2 measured from symmetric pixels, scaled to the 
shot-noise level, y is the far-field vertical position along the maximum gain circle. Gray lines: without compensation crystal, 
dashed gray line AA = lnm, solid gray line AA = 0.5nm. Black lines: with optimal compensation crystal, solid black line 
AA = lnm, dashed black line AA = 8nm. 
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narrow-band temporal and spatial filtering, and/or by using crystals that exhibit a smaller amount of walk-off. Figure 
[8] details the role of the compensation crystal. It plots the noise in the difference between Stokes parameters measured 
from symmetric pixels as a function of the vertical coordinate y along the circle corresponding to the maximum gain 
at the degenerate frequency. 



2. Broad-band frequency filtering results 
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FIG. 9: Broad-band frequency filtering. Distribution of the noise in the difference between 5*2 (5*3) measured from symmetric 
portions of the beam cross-section, scaled to the shot-noise level. In part a) a frequency filter AA = 20 nm wide, centered 
around the degenerate frequency is used; in b) AA = 60 nm. 

The results described in section llll D ll were obtained by using relatively narrow frequency filters (5-8 nm). Remark- 
ably, when a broader frequency filter is employed, the regions where Stokes parameter are correlated stretch to form 
a ring-shaped region around the pump direction (see Fig. GJJ. This kind of shape can be understood by considering 
the geometry of the downconversion cones emitted at the various frequencies by a BBO crystal. Figure ITHI is a polar 
plot of the phase matching curves (geometrical loci of the phase matched modes), with being the polar angle from 
the pump direction of propagation (z axis) and <j) the azymuthal angle around z. In this plot the same color identifies 
the same emission wavelength; dark/light thick curves correspond to two conjugate wavelengths, while the thin black 
curves are the two emission cones at the degenerate wavelength. The signal (ordinary) wave emission curves are those 
in the upper half of the plot. When considering the intersection points of a dark circle with a light circle, which 
correspond to two conjugate wavelenghts, we cannot expect any kind of entanglement, since photons arriving in these 
positions are clearly distinguishable by their different frequencies. Let us consider, instead, one of the intersection of 
the light curves (e.g the one pointed by the arrow in the plot). Here ordinary and extraordinary photons arrive with 
identical probability, and have the same wavelength. As a consequence, the photon polarization is undetermined, and 
the light is completely unpolarized. However, each time an ordinary (extraordinary) photon arrives at this position, 
an extraordinary (ordinary) photon, at the conjugate wavelength, will be found at the symmetric position. This corre- 
sponds to the intersection of the two dark curves, indicated in the plot by the second arrow. Hence, when considering 
photodetection from the two regions indicated by the arrows in the plot we can expect a high degree of polarization 
entanglement. The same reasoning can be made for any of intesection of circles corresponding to the same wavelength 
(light with light, dark with dark and thin black with thin black). By connecting all these points toghether, we can 
for example recognize the geometrical shape of the dark regions in Fig|^, where a high degree of correlation in all 
the Stokes parameters exists. By including more frequencies, the ring-shaped region if Fig. [!J> shows up. 
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FIG. 10: Polar plot of phase matching curves in a 2 mm BBO crystal cut at 49.6 degrees. 6 is the polar angle from the pump 
direction of propagation; 4> is the azymuthal angle around the pump. Black thin curves: \ S ignai = Xidler = 702nm Dark thick 
curves: \ a i gna i = Xidler = 692nra. Light thick curves:A s i gna i = Xidier = 712.29 nm. The signal is the ordinary (slow) wave, and 
correspond to curves in the upper half of the plot. Idler is the extraordinary (fast) wave. 



From the mathematical point of view, the form of the region where Stokes parameters are quantum correlated is given 
by the solution of the equation: 

A(q, fi) = A(-q, -ft) = (78) 

with x = 27r<f/(A/). By introducing the explicit form the the phase mismatch function 1241 (|78() is the equation of 
an ellipse (with a small eccentricity) centered around the origin, that is, the best correlated modes form a slightly 
asymmetric cone around the pump direction. 



IV. POLARIZATION CORRELATION: QUANTUM STATE FORMALISM 

This section is devoted to the discussion of the problem in terms of an equivalent quantum state formalism. Our 
aim is, on the one side, to give an alternative and istructive point of view on the problem, which can be compared 
with already existing quantum-state description of the problem (see e.g |lOj|). On the other side, we think that this 
section will show how the quantum field formalism developed in the first part of this paper is more powerful and 
strightforward, in terms of calculation efforts, than the commonly used quantum state formalism, at least for this 
kind of multi-mode problems. 

Equation l|ll|) defines a linear transformation acting on field operators, that maps field operators at the entrance 
face of the crystal into those at the output face. An equivalent transformation, acting on the quantum state of the 
signal/idler fields at the crystal input and mapping it into the state at the crystal output, is derived in details in 
Appendix B. As described in the appendix, in order to avoid formal difficulties coming from a continuum of modes, 
we have introduced a quantisation box both in the transverse spatial domain and in the temporal domain, so that 
the continuum of spatio-temporal modes q, fl is replaced by a discrete set of modes. 

When at the input of the parametric crystal there is the vacuum state for both signal and idler fields 

Mi„ = I"*) = II 1°; ~ n )o 10; q, n) e , (79) 

we find that the output state takes the form 

l^>e™t = II | E |n; q, fi) |n; -q, -Q) e 1 (80) 

q,n ln=0 J 

cM O) = [tanhr® ft)]" e**^ = ^f^fe^ , (81) 

coshr{q,il) \U (q,Cl)\ + 
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where the notation \n; q, f2) / e indicates the Fock state with n photons in mode (q, ft) of the ordinary /extraordinary 
polarized beam. Here the functions U ,V e are the coefficients of the operator transformation Ijlljl. and functions r 
and ip are defined by Equations I|17I18|I toghether with EH • 

The state (1811) is clearly entangled (non factorizable) with respect to the ordinary and extraordinary polarized beam 
components. 

Let us focus on two conjugate modes q, Q and — q, — fi for both the ordinary and extraordinary field components. 
These can be for example observed by using a narrow filter around the degenerate frequency f2 = and collecting 
light from two diaphragms placed around two symmetric regions in the far field zone. For brevity of notation, let us 
label these modes with the 3D vectors £ = (q x ,q y ,£l), — £ = (—q x ,—q y , — fi). When restricted to these modes, the 
state takes the form 




Cn'(-£) 



m=0 



(82) 
(83) 
(84) 



where the last two lines have been obtained by changing the dummy summation variables n, n' into m = n, N = n+n' . 
The state can be represented as a superposition of states with a fixed total number of photons N. In each N-photon 
state described by Eq. (|Sl|l 



7N,m{i) = Cm{£)CN-m(-£) 

[tanhr(0] m [tanhr(-£)] Ar - r 
coshr(£) coshr(— £) 



. e 2iJVV(-|) e 2im[V(|)-V(-D] 



(85) 



represents the probability amplitude of finding m ordinary photons, N-m extraordinary photons in mode £, and N-m 
ordinary photons, m extraordinary photons in the conjugate mode — £. The description of the state given by Eas. (|82t 
I84[) is a generalization of that derived in e.g. |ldj| . The main improvement is that our description includes the effects 
of spatial and temporal walk-off, and allows the quantitative evaluation of all the quantities of interest by using the 
parameters of a real crystal. Remarkably, when the spatial and temporal walk-off are not taken into account, it holds 
the symmetry (<f, fl) — > (—q, In this case, in Ea. (|55|) we would have r(— £) = r(£) and £) = and all the 

coefficients 7Ar jm (£) would be identical for a given N, so that all the terms in the expansion would have the same 
weight, thus leading to a "maximally entangled state for notarization? [TfT| . 

Coming to Stokes parameter correlation we notice the following property of the state: 



N 



At&M?) - 4(lMe(l)] \^)% = E c m (|)c^_ m (-0(2m - N) 



m=0 



N -m;£ 



N 



-|)(86) 
(87) 



By recalling the definition of the Stokes operator densities given by Eas. l|38l4()|) ai(£) = Al(£)A (£) — Al(£)A e (£), 
0-2(1) = Al(£)A e (£) + Al(£)A (£) and <7 3 (|) = -«[ij(|)i e (|) - A|(|)i (|)], we can hence conclude that the state is 
an eigenstate of <7i(£) + &x{—£) with zero eigenvalue. On the other side, we have 



Al(£)Ae& 



N-l 



c m{£) c N-m(-£W(™ + 1)CW - m) 



m— 



m + 1 ; £ 



N-m-l;£ 



N - m; -£ 



to; — £ 



m— 1 
N-l 



= c I+ iK> w _,-i(-ovo + W - o 



to;£ 



N -m;£ 



iV — m + 1; — £ 



N-l-l;Z 



1=0 



m-l;-^9) 
(90) 
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where the last line has been obtained by introducing the summation index I = m — 1. This implies that the equation 

4(1)^(1) 

is verified if and only if 

( -Hi 

for all N = 0, +oo and m = 0,7V - 1. 
Similar considerations for the hermitian conjugate operator A\{S)A (^) — A\{— ^)A {— £) lead to the equivalent con- 
dition 

Cm(£,)cN-m{-£,) = C m _l(|)cjV_ m +l(-|) (93) 

for all iV = 0, +00 and m = 1, TV. 

Hence, the state is also an eigenstate of both (72(C) — 02 ( — £) and <73(£) — 03(— £), with zero eingevalue, if and only if 
the conditions (|92I93I ) are satisfied. These conditions amount to requiring that all the coefficients in the expansion of 
the N-photon state 1)84(1 are identical, and that the N photon state is a superposition with equal probability amplitude 
of all the possibile partitions in m ordinary and N-m extraordinary photons (m=0,N) in mode £, with N-m ordinary 
and m extraordinary photons in the conjugate mode — £. This is the mathematical equivalent of the commonly used 
statement 11 ordinary and extraordinary photons in mode £ are not distinguishable, but each time we have m ordinary 
and N-m extraordinary photon in mode £, there are N-m ordinary and m extraordinary photons in mode — £" . For 
modes having a non vanishing parametric gain the conditions (I92I93[) amount to requiring 

tanhr(f)e 2i ^ = tanhr(-f)e 2iV,( -^ , (94) 

a condition that is satisfied only in the presence of the symmetry A(q, Q) = A(— q, — fi). This in turns implies the 
absence of spatial walk-off between the two waves(i.e. the two modes correspond to the intersection of the down- 
conversion cones) and the absence of temporal walk-off (use of a narrow frequency filter and/or compensation by 
means of a second crystal). 
Formula can be also written as: 

U$V*(-$ = U(-0V*(0 . (95) 

By comparing with equation (|73|l . we notice that this is the condition that ensures that the correlation between 
Stokes parameter measured from symmetric pixels calculated in Section IlII Dl reaches its maximum value. Hence, in 
the framework of the quantum state formalism, we start to recover the same results of Section IlII Bl as it obviously 
must be. One could proceed further on, and derive quantitative results for the correlation, as those showed by Figs l4l9l 
but at this point it should be rather clear (and for sure we are not the first ones to notice this) how the quantum state 
formalism, although instructive, is cumbersome and not transparent in comparison with the quantum field formalism. 

V. CONCLUSIONS 

In conclusion, we have shown that the polarization entanglement of photon pairs emitted in parametric down- 
conversion survives in high gain regimes, where the number of converted photons can be rather large. In this case, it 
takes the form of non-classical spatial correlations of all light Stokes operators associated to polarization degrees of 
freedom. We have shown that in the regions where the two rings intersect (in a ring-shaped region around the pump 
direction when a broad frequency filter is employed) all the Stokes operators are highly correlated at a quantum level, 
realizing in this way a macroscopic polarization entanglement. Although Stokes parameters are extremely noisy and 
the state is unpolarized, measurement of a Stokes parameter in any polarization basis in one far-field region determines 
the Stokes parameter collected from the symmetric region, within an uncertainty much below the standard quantum 
limit. 

We call this situation "polarization entanglement" because, on the one side, the quantum state derived in Section llVl is 
entangled with respect to polarization degrees of freedom, and, on the other, because in our description there is no gap 
in the passage from the single photon pair regime, where the polarization entanglement is a widely accepted concept, 
to the multiple photon pair regime. However, we want to remark that for spontaneous parametric down-conversion 
there is no way, to our knowledge, to derive a sufficient criterium for inseparability based on the degree of correlation 



(91) 



-£) = Cm+l(O c iV-m-l(-|j 



(92) 
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of the Stokes operators, as this derived in [26j and generalized by |3[. This depends on the fact that the average 
values of commutators (and anticommutators) of Stokes operators are in this system intrinsically state dependent, 
at difference to what happens in the experiment performed in 0, where bright entangled beams were used. Further 
discussion about this important point is postponed to a future publication. 

We have developed a multi-mode model for spontaneous parametric down-conversion, both within the framework 
of quantum field formalism and quantum state formalism. They are valid in any gain regime, from the single photon 
pair production to the high gain regime where the number of downconverted photons can be rather large. The model 
allows quantitative estimations of all the quantities of interest, by using empirical parameters of real crystals. We 
hope that this description can be a useful tool for experimentalists working in this field. 

Quite interesting, and to our knowledge completely novel, are the results concerning the correlation of Stokes 
parameters observed by using a broad frequency filter, described in Section IIII D 21 They basically show how by 
increasing the number of temporal degrees of freedom in play , the number of spatial degrees of freedom which are 
simultaneously entangled increases, so that the two isolated correlated spots in Figure 0] become the ring shaped 
region of Figure where many symmetric spots are correlated in pairs. 

We believe that this form of entanglement, with its increased complexity in terms of degrees of freedom (photon 
number, polarization, temporal and spatial degrees of freedom)can be quite promising for new quantum information 
schemes. 



Appendix A 



In this Appendix we calculate the phase shift induced by the propagation of the down-converted fields through a 
compensation crystal, and we evaluate the length of this second crystal necessary for optimal walk-off compensation. 
As shown by the scheme of FigEI we assume that after producing down-conversion in a first crystal (BBOl), the 
pump beam is eliminated. The polarisations of the downconverted beams is then rotated by 90 degrees, and they 
pass through a second crystal(BB02) of length l' c , identical to the first one. 

In the region between the second crystal and the lens L the ordinary/extraordinary field operators can be written as: 



A (q,n,z) = i™ t (g,n)exp[ifc 02 (g,n)^]exp[i^ ac ^-^)] 
A e (q,n,z) = i° ui (g,«)exp[ifc ez (g,fi)^]ex P [i^ ac (z-0] 



(96) 
(97) 

The first phase shift accounts for propagation inside the compensation crystal. Here k oz (q,Q), k ez (q,Q) are the 
projections along z-axis of the ordinary/extraordinary wave-vectors inside the crystal, whose explicit expressions 
depend on the linear properties of the crystal as described by Eq.J33J). The second phase shift accounts for paraxial 

2 

propagation in vacuum 4> vac {z) = (k — %c)z , fc = 27r/A 



In the far field plane, all the results described in Sections IIII CI ITTl Dl remain unchanged provided that one makes the 
following substitutions: 

(98) 
(99) 
(100) 
(101) 

where global phase factors have been omitted, since they do not affect the results. 

This transformation leaves unchanged all the results described in Sec. IIII CI (noise and correlation for measurements 
of Stokes operators and 1). For the second and third Stokes parameters (Sec. IIII D|l , while the transformation 
does not affect the amount of noise of the measurement, given by Eqs. 1691711 it does affect the correlation between 
measurements from symmetric pixels (Eqs. I70I72|I 

^2Re{w*(-f,-w)W(f,w)V*(-f,-tj)V(f,w)e i ^ (:? '^} , (102) 



U (x,il) - 


+ U e (q, fi)exp 




1= 3 TJ 


V (x,Q) - 


+ V e (q, fi)exp 






u e (x,n) - 


+ U (q,£l)exp 


[ik ez (q^)Q 




v (£,n) - 


+ V (q, fi)exp 







H 2 {x, 0) 

cj> c {x,Vi) = [k ez (q,Q) + k oz (-q, -O) - k oz (q,Q) - k ez (-q, -to)]l' c 

t 



[A(-q,-Q)-A(q,n)}f 



q=x%f 



(103) 
(104) 
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On the other side, by using the explicit expression of the gain functions in Eqs. (|19I20|1 . we have 

arg {U* (-x, -lj)U{x, uj)V* {-x, ~lu)V(x, u)} = 2tp(q, O) - 2^(-q, -CI) (105) 

q=x^ 

with 

W,0) = taa-^ani ^ff } (106) 
tanhfT 

w A(£fi)-^— . (107) 

The last line has been obtained by taking the limit A(g, fi) <C 1; this is meaningful since the most important 
contribution to the correlation function is given by phase matched modes. The phase factor (|105[) can be partially 
compensated by the phase shift induced by propagation in the second crystal l|104f) . Best compensation is achieved 
for 

I' tanhcr 

i = — ^ 

In this conditions the value of the correlation between measurements from symmetric pixel (the value of the function 
H2) is maximized by the presence of a compensation crystal. 



Appendix B 



Equation l|ll|) defines a linear transformation acting on field operators, that maps field operators at the entrance 
face of the crystal into those at the output face. The aim of this appendix is to find is an equivalent transformation 
acting on the quantum state of the signal/idler fields at the crystal input and mapping it into the state at the crystal 
output. 

In order to avoid formal difficulties coming from a continuum of modes, we introduce a quantisation box of side b in 
the transverse plane, with periodic boundary conditions. In this way the continuum of wave- vectors (f is replaced by a 
set of discrete wave vectors = (l x u x + lyU y )^- l x ,l y = 0,±1,±2.... In the same way, we introduce a quantization box 
in the time domain of length T, with periodic boundaries, so that we need to consider only a discrete set of temporal 
frequencies fl p = p^- p — 0, ±1... The free field commutation relation (|12[) are thus replaced by their discrete version 

Clp),Aj(q^,Cl s ) = i5(,j^,mji„m»i5n,s ,i,3 = 0,e. 
For brevity of notation, in the following we shall indicate the spatio-temporal mode qk Q p with the three dimensional 

vector £, and we shall not write explicitly the modal indices. 

The input/output transformation 111(1 can be written in a equivalent way as: 



with 



and 



R — R0R1R2 , 



Ro — 



Ri = 



R-2 



exp I $>(D \AKoAK-o - i (i)i e (-D" 



exp \ ij^m [Aims + aksM-z) 



(109) 
(110) 

(111) 
(112) 
(113) 
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with functions ¥>(|),r(|), 0(|) defined h V Ec l s - D2d> toghether with (PEOIEJ- 

In order to demonstrate the ansatz l|109fl . we first notice that the action of operators Ro and R2 on field operators 
corresponds to phase rotations. For any operator c, for which [c, c^] = 1, we have 

e -is C t Cce i sc + c = e is c _ ^ 114 j 

As a consequence, 

i?Si o (£)i? = i (|y[^>^] (115) 

RlA e (£)R = Aeffijl^-d-rt-d] . (116) 

Operator R\ is the product of an infinity of two mode squeezing operators, each of them acting on the couple of 
modes (£) in the signal beam and (— £) in the idler beam. For any couple of independent boson operators c%, C2, and 
for r real, it holds the identity 

Hence, letting ci — > A D (£), C2 — ► A e (— £), we have 

i?Ii (e)i?i = A (0 cosh r(|) + A\ (-f)sinhr(e) (118) 
flfi e (C)4 = Mi) coshr(-f) + it(-e) sinhr(-e) . (119) 

Finally, by letting also the operator R2 act: 

£Ji?}i?Ji o (0i? 44 - e'^^+^^l {i o (0coshr(e> ie( « ) +it(-Osinhr(e)e~ ie( « ) } (120) 
= + At(-t)V(0} (121) 

where in passing from the first to the second line we used the relation l|18|) , which is a consequence of the unitarity of 
the transformation Qllfl. Moreover, we have 

RlRlRlA.iORQR^ = jlMS-rt-d] {i e (|)coshr(-|)e ie ^ 

+ it(-|)sinhr(-De- ie (-^} (122) 

= e -w-s {ie(e>(-l) + i f (-IM-D} • (123) 

Finally, taking into account the relation 117|) . which is again a consequence of the unitarity of the transformation i|ll[) . 
we recover the input/output transformation 1111) . 

Any quantum mechanical expectation value of the output operators (mean values, correlation functions etc.) taken 
on the input state, is equivalent to the quantum mechanical expectation value of the input operators taken on the 
transformed state: 

Wont = R Win ( 124 ) 

In the following we shall derive the form of the output state, when at the input of the parametric crystal there is the 
vacuum state for both signal and idler fields. 

l^) in = l^c)=Jj|0;l) o |0;-|) e , (125) 
f 

where the notation n; O indicates the Fock state with n photons in mode (£) of the ordinary/extraordinary 

/ o/e 

polarized beam. 

First of all we notice that the operator Ri has no effect on the vacuum state, corresponding to a phase rotation of 
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the vacuum. For what concern operator Hi, by using proper operator ordering techniques (see e.g. |25|) pag. 75), it 
can be recasted in the following form (disentangling theorem) 



Ri = I \ < < 

G(0 = tanh[r(|)] 
g(0 = log{cosh[r(0]} 

By letting this operator acting on the vacuum state 



jjj e G(f)At(ai!(-a e -s(a^^ 



Ri \vae) = I] — TTTTfT E [ tan M0 
- cosh[r(OJ n=0 L 



n 








n;|) 


n; 









(126) 

(127) 
(128) 



(129) 



where the usual expansion of the exponential operator, expM — X^^Lo ITF' nas been used , toghether with the 
standard action of boson creation operators on Fock states. Finally, by adding the action of operator Rq, 



R Rt \vac) = [] ) Y, ftanhr(C) 
- cosh[r(OJ n L 



n;-£ 



(130) 



the output state can be written in the form 



Wont = n{E C ««) 



c„(0 = 



1 

coshr(£) 



tanhr(£) 



3 2inV(?) _ 



^o(l)Ve(-|) 



Uo(0 



2n+l 



(131) 
(132) 
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